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7. Thermoacoustics 


Thermodynamic and fluid-dynamic processes in 
sound waves in gases can convert energy from one 
form to another. In these thermoacoustic processes 
[7.1,2], high-temperature heat or chemical energy 
can be partially converted to acoustic power, 
acoustic power can produce heat, acoustic power 
can pump heat from a low temperature or to a high 
temperature, and acoustic power can be partially 
converted to chemical potential in the separation 
of gas mixtures. In some cases, the thermoacoustic 
perspective brings new insights to decades-old 
technologies. Well-engineered thermoacoustic 
devices using extremely intense sound approach 
the power conversion per unit volume and the 
efficiency of mature energy-conversion equipment 
such as internal combustion engines, and the 
simplicity of few or no moving parts drives the 
development of practical applications. 

This chapter surveys thermoacoustic energy 
conversion, so the reader can understand how 
thermoacoustic devices work and can estimate 
some relevant numbers. After a brief history, an 
initial section defines vocabulary and establishes 
preliminary concepts, and subsequent sections 
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explain engines, dissipation, refrigeration, and 
mixture separation. Combustion thermoacoustics 
is mentioned only briefly. Transduction and 
measurement systems that use heat-generated 
surface and bulk acoustic waves in solids are not 
discussed. 



7.1 History 

The history of thermoacoustic energy conversion has 
many interwoven roots, branches, and trunks. It is 
a complicated history because invention and technology 
development outside of the discipline of acoustics 
have sometimes preceded fundamental understand¬ 
ing; at other times, fundamental science has come 
first. 

Rott [7.3, 4] developed the mathematics describ¬ 
ing acoustic oscillations in a gas in a channel with an 
axial temperature gradient, with lateral channel dimen¬ 
sions of the order of the gas thermal penetration depth 
(typically ~ 1 mm), this being much shorter than the 
wavelength (typically ~lm). The problem had been 
investigated by Rayleigh and by Kirchhoff, but without 
quantitative success. In Rott’s time, the motivation to un¬ 


derstand the problem arose largely from the cryogenic 
phenomenon known as Taconis oscillations - when 
a gas-filled tube is cooled from ambient temperature 
to a cryogenic temperature, the gas sometimes oscillates 
spontaneously, with large heat transport from ambient to 
the cryogenic environment. Yazaki [7.5] demonstrated 
convincingly that Rott’s analysis of Taconis oscillations 
was quantitatively accurate. 

A century earlier, Rayleigh [7.6] understood the 
qualitative features of such heat-driven oscillations: “If 
heat be given to the air at the moment of greatest con¬ 
densation (i.e., greatest density) or be taken from it at 
the moment of greatest rarefaction, the vibration is en¬ 
couraged.” He had studied Sondhauss oscillations [7.7], 
the glassblowers’ precursor to Taconis oscillations. 


Part B 




















Part BI 7.2 


240 Part B 


Physical and Nonlinear Acoustics 


Applying Rott’s mathematics to a situation where the 
temperature gradient along the channel was too weak to 
satisfy Rayleigh’s criterion for spontaneous oscillations, 
Hofler [7.8] invented a standing-wave thermoacoustic 
refrigerator, and demonstrated [7.9] again that Rott’s 
approach to acoustics in small channels was quan¬ 
titatively accurate. In this type of refrigerator, the 
coupled oscillations of gas motion, temperature, and 
heat transfer in the sound wave are phased in time so 
that heat is absorbed from a load at a low tempera¬ 
ture and waste heat is rejected to a sink at a higher 
temperature. 

Meanwhile, completely independently, pulse-tube 
refrigeration was becoming the most actively investi¬ 
gated area of cryogenic refrigeration. This development 
began with Gifford’s [7.10] accidental discovery and 
subsequent investigation of the cooling associated with 
square-wave pulses of pressure applied to one end of 
a pipe that was closed at the other end. Although 
the relationship was not recognized at the time, this 
phenomenon shared much physics with Holler’s refrig¬ 
erator. Mikulin ’s [7.11] attempt at improvement in heat 
transfer in one part of this basic pulse-tube refrigerator 
led unexpectedly to a dramatic improvement of perfor¬ 
mance, and Radebaugh [7.12] realized that the resulting 


orifice pulse-tube refrigerator was in fact a variant of the 
Stirling cryocooler. 

Development of Stirling engines and refrigerators 
started in the 19th century, the engines at first as an 
alternative to steam engines [7.13]. Crankshafts, multi¬ 
ple pistons, and other moving parts seemed at first to be 
essential. An important modern chapter in their develop¬ 
ment began in the 1970s with the invention of free-piston 
Stirling engines and refrigerators, in which each pis¬ 
ton’s motion is determined by interactions between the 
piston’s dynamics and the gas’s dynamics rather than 
by a crankshaft and connecting rod. Analysis of such 
complex, coupled phenomena is complicated, because 
the oscillating motion causes oscillating pressure dif¬ 
ferences while simultaneously the oscillating pressure 
differences cause oscillating motion. Ceperley [7.14,15] 
added an explicitly acoustic perspective to Stirling en¬ 
gines and refrigerators when he realized that the time 
phasing between pressure and motion oscillations in 
the heart of their regenerators is that of a traveling 
acoustic wave. Many years later, acoustic versions of 
such engines were demonstrated by Yazaki [7.16], de- 
Blok [7.17], and Backhaus [7.18], the last achieving 
a heat-to-acoustic energy efficiency comparable to that 
of other mature energy-conversion technologies. 
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7.2.1 Pressure and Velocity 

For a monofrequency wave, oscillating variables can be 
represented with complex notation, such as 

p(x,t) = p m + Re[pi(x)e lcor ] (7.1) 


as two first-order differential equations coupling two 
variables, pressure p\ and the x component of volume 
velocity, U\\ 


dpi 

dx 


\COPxxx 

A 



(7.3) 


for the pressure p , where p m is the mean pressure, Re (z) 
indicates the real part of z, co = In f is the angular fre¬ 
quency, / is the frequency, and the complex number p \ 
specifies both the amplitude and the time phase of the 
oscillating part of the pressure. For propagation in the x 
direction through a cross-sectional area A in a duct, p\ is 
a function of x. In the simple lossless, uniform-area situ¬ 
ation the sinusoidal x dependence can be found from the 
wave equation, which can be written with ico substituted 
for time derivatives as 


2 , 2 d PI 

" p ‘ +c 



(7.2) 


where c 2 = (dp/ 3p) s is the square of the adiabatic sound 
speed, with p the density and s the entropy. In thermo¬ 
acoustics, intuition is well served by thinking of (7.2) 


dU\ _ icoA 
dx PmC 2 ^ 1 


(7.4) 


For a reasonably short length of duct Ax, these can be 
written approximately as 


Api = —i coL U 1 , 


P 1 = - 


1 

i^C 


A Ui , 


where 


L = 
C = 


Pm Ax 

A 

A Ax 
PmC 2 


(7.5) 

(7.6) 


(7.7) 

(7.8) 
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Table 7.1 The acoustic-electric analog. The basic building 
blocks of thermoacoustic energy-conversion devices are in- 
ertances L, compliances C, viscous and thermal resistances 
R y isc and R± eTm , gains (or losses) G, and power transducers 


Acoustic variables 

Electric variables 

Pressure p\ 

Voltage 

Volume velocity U\ 

Current 

Inertance L 

Series inductance 

Viscous resistance R v i sc 

Series resistance 

Compliance C 

Capacitance to ground 

Thermal-hysteresis 
resistance R t herm 

Resistance to ground 

Gain/loss along 

Proportionally controlled 

temperature gradient, G 

current injector 

Stroke-controlled transducer 

Current source 

Force-controlled transducer 

Voltage source 


As introduced in Chap. 6 (Physical Acoustics) (7.3) is 
the acoustic Euler equation, the inviscid form of the mo¬ 
mentum equation. It describes the force that the pressure 
gradient must exert to overcome the gas’s inertia. In the 
lumped-element form (7.5), the geometrical and gas- 
property aspects of the inertia of the gas in the duct are 
represented by the duct’s inertance L. Similarly intro¬ 
duced in Chap. 6, (7.4) and (7.6) are, respectively, the 
differential and lumped-element versions of the continu¬ 
ity equation combined with the gas’s adiabatic equation 
of state. These describe the compressibility of the gas 
in response to being squeezed along the x-direction. In 
the lumped-element form (7.6), the geometrical and gas- 
property aspects of the compressibility of the gas in the 
duct are represented by the duct’s compliance C. 

Accurate calculation of the wave in a long duct or 
a long series of ducts requires analytical or numerical 
integration of differential equations [7.1], but for the pur¬ 
poses of this chapter we will consider lumped-element 
approximations, because they are satisfactorily accurate 
for many estimates, they allow a circuit model repre¬ 
sentation of devices, and they facilitate intuition based 
on experience with alternating-current (AC) electrical 
circuits. Table 7.1 shows the analogy between acoustic 
devices and electrical circuits. 

For example, the closed-open resonator shown in 
Fig. 7.1a has a resonance frequency that can be calcu¬ 
lated by setting its length Ai reso equal to a quarter 
wavelength of sound; the result is / res0 = c/4Ai reso . The 
simplest lumped-element approximation of the quarter- 
wavelength resonator is shown in Fig. 7.1b. We assign 
compliance C reS o to the left half of the resonator, be- 
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b) 
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Fig.7.1a-d Quarter-wavelength resonator, (a) A quarter- 
wavelength resonator: a tube closed at one end and open 
at the other. (b) A simple lossless lumped-element model 
of the resonator as a compliance C reso in series with an 
inertance L reso . (c) A more complicated lumped-element 
model, including thermal-hysteresis resistance Rthevm in the 
compliance, viscous resistance R v i sc in the inertance, and 
radiation impedance /? ra d + icoL m &. (d) An illustration of the 
first law of thermodynamics for a control volume, shown 
enclosed by the dashed line, which intersects the stack of 
a well-insulated standing-wave engine. In steady state, the 
thermal power Qn that flows into the system at the hot heat 
exchanger must equal the total power H that flows along 
the stack 


cause the compressibility of the gas is more important 
than its inertia in the left half, where \p\ | is high and 
| Gi| is low. Similarly, we assign inertance L reso to the 
right half of the resonator, because the inertia is more 
important in the right half where \U\\ is high. Setting 
Ax = Ax r eso/2 in (7.7) and (7.8), and recalling that 
the resonance frequency of an electrical LC circuit is 
given by or = 1/LC, we find / res0 = c/7rAx reS o> dif¬ 
fering from the exact result only by a factor of 4 /tt. 
Such accuracy will be acceptable for estimation and 
understanding in this chapter. 
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Circuit models can include as much detail as is nec¬ 
essary to convey the essential physics [7.1]. When the 
viscosity is included in the momentum equation, the 
lumped-element form becomes 

A/?i = - (i coL + R visc ) U\ , (7.9) 

and when the thermal conductivity and a temperature 
gradient in the x direction are included in the continuity 
equation, the lumped-element form becomes 

AU\ = — (icoC + —- j p\ + GT/i n? i . (7.10) 

Figure 7.1c is a better model than Fig. 7.1b for the 
closed-open resonator of Fig. 7.1a, because it includes 
thermal-hysteresis resistance in the compliance, viscous 
resistance in the inertance, and the inertial and resistive 
radiation impedance at the open end. This model would 
yield reasonably accurate estimates for the free-decay 
time and quality factor of the resonator as well as its 
resonance frequency. 

Ducts filled with porous media and having tempera¬ 
ture gradients in the x-direction play a central role in 
thermoacoustic energy conversion [7.1]. The pore size 
is characterized by the hydraulic radius rh, defined as the 
ratio of gas volume to gas-solid contact area. In a circu¬ 
lar pore, rh is half the circle’s radius; in the gap between 
parallel plates, rh is half the gap. The term stack is used 


for a porous medium whose rh is comparable to the gas 
thermal penetration depth 

^therm — -yj 2k/ Ojp m Cp , (7.11) 

where k is the gas thermal conductivity and c p is its 
isobaric heat capacity per unit mass, while a porous 
medium with rh <<C <$therm is known as a regenerator. 
The viscous penetration depth 

^visc — y/ 2/x /CO Pm ? (7.12) 

where /z is the viscosity, is typically comparable to, 
but smaller than, the thermal penetration depth. If the 
distance between a particular mass of gas and the nearest 
solid wall is much greater than the penetration depths, 
thermal and viscous interactions with the wall do not 
affect that gas. 

In such porous media with axial temperature gradi¬ 
ents, the gain/loss term GU m ,\ in (7.10) is responsible 
for the creation of acoustic power by engines and for the 
thermodynamically required consumption of acoustic 
power by refrigerators. This term represents cyclic ther¬ 
mal expansion and contraction of the gas as it moves 
along and experiences the temperature gradient in the 
solid. In regenerators, the gain G is nearly real, and 
AU\ caused by the motion along the temperature gradi¬ 
ent is in phase with U\ itself, because of the excellent 
thermal contact between the gas in small pores and the 


Table 7.2 Expressions for the lumped-element building blocks L, R v i sc , C, 7? t herm, an d GU i, and for the total power H, 
in the boundary-layer limit rh ^ 8 and in the small-pore limit r^<^ 8. The symbol in the small-pore limit indicates 
that the numerical prefactor depends on the shape of the pore 
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walls of the pores. (Positive G indicates gain; negative 
G indicates loss.) In stacks, a nonzero imaginary part 
of G reflects imperfect thermal contact in the pores and 
the resultant time delay between the gas’s cyclic mo¬ 
tion along the solid’s temperature gradient and the gas’s 
thermal expansion and contraction. 

Table 7.2 gives expressions for the impedances in 
the boundary-layer limit, r h ^ <$ t herm an d Hi ^ <$visc , 
which is appropriate for large ducts, and in the small- 
pore limit appropriate for regenerators. Boundary-layer- 
limit entries are usefully accurate in stacks, in which 
rh ~ 8. General expressions for arbitrary r h and many 
common pore shapes are given in [7.1]. 

The lumped-element approach summarized in Ta¬ 
ble 7.1 and the limiting expressions given in Table 7.2 
are sufficient for most of this overview, but quantitatively 
accurate thermoacoustic analysis requires slightly more 
sophisticated techniques and includes more phenom¬ 
ena [7.1]. Every differential length dx of duct has dL, 
dC, d^ v i sc , and d(l/R t hermX and if d T m /dx 7^ 0 it also 
has dG, so a finite-length element is more analogous 
to an electrical transmission line than to a few lumped 
impedances. In addition to smoothly varying x depen¬ 
dencies for all variables, important phenomena include 
turbulence, which increases R w i sc above the values given 
in Table 7.2; pore sizes which are in neither of the limits 
given in Table 7.2; nonlinear terms in the momentum and 
continuity equations, which cause frequency doubling, 
tripling, etc., so that the steady-state wave is a superposi¬ 
tion of waves of many frequencies; and streaming flows 
caused by the wave itself. Many of these subjects are 
introduced in Chap. 8 (Nonlinear Acoustics). Thermo¬ 
acoustics software that includes most or all of these 
phenomena and has the properties of several commonly 
used gases is available [7.19,20]. 

For estimating the behavior of thermoacoustic de¬ 
vices, it is useful to remember some properties of 
common ideal gases [7.21]. The equation of state is 


P = 


pR umv T 

M 


(7.13) 


where R un iv = 8.3 J/(moleK) is the universal gas con¬ 
stant and M is the molar mass. The ratio of isobaric to 
isochoric specific heats, y, is 5/3 for monatomic gases 
such as helium and 7/5 for diatomic gases such as ni¬ 
trogen and air near ambient temperature, and appears in 
both the adiabatic sound speed 


c = 


yRunivT 


M 


(7.14) 


and the isobaric heat capacity per unit mass 



yR 

univ 

(y — M ' 


(7.15) 


The viscosity of many common gases (e.g., air, helium, 
and argon) is about 


/ t \ °-7 

M-px 10~ 5 kg/ms) (^qqj^ ) ’ (7 - 16) 

and the thermal conductivity k can be estimated by 
remembering that the Prandtl number 


/ZCp 

k 


(7.17) 


is about 2/3 for pure gases, but somewhat lower for gas 
mixtures [7.22]. 


7.2.2 Power 


In addition to ordinary acoustic power E , the time- 

• • 

averaged thermal power Q, total power H , and exergy 
flux X are important in thermoacoustic energy con¬ 
version. These thermoacoustic powers are related to 
the simpler concepts of work, heat, enthalpy, and ex¬ 
ergy that are encountered in introductory [7.21] and 
advanced [7.23] thermodynamics. 

Just as acoustic intensity is the time-averaged prod- 
uct of pressure and velocity, acoustic power E is the 
time-averaged product of pressure and volume velocity. 
In complex notation, 

E = ^Re (p\U\) , (7.18) 

where the tilde denotes complex conjugation. At trans¬ 
ducers, it is apparent that acoustic power is closely 
related to thermodynamic work, because a moving pis¬ 
ton working against gas in an open space with volume V 
transforms mechanical power to acoustic power (or vice 
versa) at a rate f § p dV, which is equal to (7.18) for 
sinusoidal pressure and motion. Resistances R dissipate 
acoustic power; the gain/loss term GU 1 in components 
with temperature gradients can either consume or pro¬ 
duce acoustic power, and inertances L and compliances 
C neither dissipate nor produce acoustic power, but 
simply pass it along while changing p\ or U\. 

Time-averaged thermal power Q (i.e., time-averaged 
heat per unit time) is added to or removed from the 
gas at heat exchangers, which are typically arrays of 
tubes, high-conductivity fins, or both, spanning a duct. 
Thermal power can be supplied to an engine by high- 
temperature combustion products flowing through such 
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tubes or by circulation of a high-temperature heat- 
transfer fluid through such tubes and a source of heat 
elsewhere. Thermal power is often removed from en¬ 
gines and refrigerators by ambient-temperature water 
flowing through such tubes. 

Of greatest importance is the total time-averaged 
power 



1 

2 


p m Re ( h\u{) —k 


d T m 
dx 


d A , 


(7.19) 


based on the v component u of velocity and the enthalpy 
h per unit mass, which is the energy of most utility in 
fluid mechanics. Fig. 7. Id uses a simple standing-wave 
engine (discussed more fully below) to illustrate the 
centrality of H to the first law of thermodynamics in 
thermoacoustics. The figure shows a heat exchanger and 
stack in a quarter-wavelength resonator. When heat is 
applied to the hot heat exchanger, the resonance is driven 
by processes (described below) in the stack. The dashed 
line in Fig. 7. Id encloses a volume whose energy must 
be constant when the engine is running in steady state. If 
the side walls of the engine are well insulated and rigid 
within that volume, then the only energy flows per unit 


time into and out of the volume are Qw and whatever 

power flows to the right through the stack. We define 

the total power flow through the stack to be H , and the 

• • 

first law of thermodynamics ensures that H = Qn in this 
simple engine. 

As shown in (7.19), the total power H is the sum 
of ordinary steady-state heat conduction (most impor¬ 
tantly in the solid parts of the stack or regenerator and 
the surrounding duct walls) and the total power carried 
convectively by the gas itself. Analysis of the gas con¬ 
tribution requires spatial and temporal averaging of the 
enthalpy transport in the gas [7.4], and shows that the 
most important contributions are acoustic power flowing 
through the pores of the stack and a shuttling of energy 
by the gas that occurs because entropy oscillations in the 
gas are nonzero and have a component in phase with ve¬ 
locity. Remarkably, these two phenomena nearly cancel 
in the small pores of a regenerator. 

Finally, the exergy flux X represents the rate at which 
thermodynamic work can be done, in principle, with 
unrestricted access to a thermal bath at ambient tempera¬ 
ture [7.1,23]. Exergy flux is sometimes used in complex 
systems to analyze sources of inefficiency according to 
location or process. 


7.3 Engines 

Implicit in Rayleigh’s criterion [7.6] for spontaneous 
oscillations, “If heat be given to the air at the moment 
of greatest (density) or be taken from it at the moment 
of greatest rarefaction, the vibration is encouraged,” is 
the existence of a vibration in need of encouragement, 
typically a resonance with a high quality factor. To¬ 
day, we would express Rayleigh’s criterion in either 
of two ways, depending on whether we adopt a La- 
grangian perspective, focusing on discrete masses of 
gas as they move, or an Eulerian perspective, focus¬ 
ing on fixed locations in space as the gas moves past 
them. In the Lagrangian perspective, some of the gas in 
a thermoacoustic engine must experience § p dV > 0, 
where V is the volume of an identified mass of gas. In 
the Eulerian perspective, part of the device must have 
dE/dx > 0, arising from Re ( p\ dU\/ dx) > 0 and the 
G term in (7.10). 

By engine we mean a prime mover, i.e., something 
that converts heat to work or acoustic power. We describe 
three varieties: standing-wave engines, traveling-wave 
engines, and pulse combustors. 


7.3.1 Standing-Wave Engines 

Continuing the quarter-wavelength example introduced 
in Fig. 7.1, Fig. 7.2 shows a simple standing-wave 
engine, of a type that is available as a classroom demon¬ 
stration [7.24]. A stack and hot heat exchanger are near 
the left end of the resonator, and its right end is open to 
the air. When heat Qh is first applied to the hot heat ex¬ 
changer, the heat exchanger and the adjacent end of the 
stack warm up, establishing an axial temperature gradi¬ 
ent from hot to ambient in the stack. When the gradient 
is steep enough (as explained below), the acoustic oscil¬ 
lations start spontaneously, and grow in intensity as more 

heat is added and a steady state is approached. In the 

• • 

steady state, total power H equal to Qh (minus any heat 
leak to the room) flows through the stack, creating acous¬ 
tic power, some of which is dissipated elsewhere in the 
resonator and some of which is radiated into the room. 
To the right of the stack, where an ambient-temperature 
heat exchanger would often be located, the heat is carried 
rightward and out of this resonator by streaming-driven 
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convection in the resonator, while fresh air at ambient 
temperature 7 a streams inwards. 

The most important process in the standing-wave 
engine is illustrated in Fig. 7. 2d and Fig. 7.2e from 
a Lagrangian perspective. Figure 7. 2d shows a greatly 
magnified view of a mass of gas inside one pore of 
the stack. The sinusoidal thermodynamic cycle of that 
mass of gas in pressure p and volume V is shown in 
Fig. 7.2e; the mass’s temperature, entropy, density, and 
other properties also vary sinusoidally in time. How- 


Fig. 7.2a-e A standing-wave engine, (a) A hot heat ex¬ 
changer and a stack in a quarter-wavelength resonator. Heat 
Qu is injected at the hot heat exchanger, and the device ra¬ 
diates acoustic power into the surroundings, (b) Total 
power flow H and acoustic power E as functions of position 
x in the device. Positive power flows in the positive-x direc¬ 
tion. (c) Lumped-element model of the device, (d) Close-up 
view of part of one pore in the stack of (a), showing a small 
mass of gas going through one full cycle, imagined as four 
discrete steps. Thin arrows represent motion, and wide, 
open arrows represent heat flow, (e) Plot showing how the 
pressure p and volume V of that small mass of gas evolve 
with time in a clockwise elliptical trajectory. Tick marks 
show approximate boundaries between the four numbered 
steps of the cycle shown in (d) 

ever, for qualitative understanding of the processes, we 
describe the cycle as if it were a time series of four dis¬ 
crete steps, numbered 1-4 in Fig. 7. 2d and Fig. 7.2e. 
In step 1, the gas is simultaneously compressed to 
a smaller volume and moved leftward by the wave a dis¬ 
tance 2 |§i | , which is much smaller than the length Ax 
of the stack. While the gas is moving leftwards, the 
pressure changes by 2\p\ \ and the temperature would 
change by 2 | T\ \ = 2 T m (y — 1) \p\ \ /yp m if the process 
were adiabatic. This suggests the definition of a critical 
temperature gradient 

VT’crit = |7i| / |^i| . (7.20) 

Thermal contact in the stack’s large pores is imperfect, 
so step 1 is actually not too far from adiabatic. However, 
the temperature gradient imposed on the stack is greater 
than VT cr it, so the gas arrives at its new location less 
warm than the adjacent pore walls. Thus, in step 2, heat 
flows from the solid into the gas, warming the gas and 
causing thermal expansion of the gas to a larger volume. 
In step 3, the gas is simultaneously expanded to a larger 
volume and moved rightward by the wave. It arrives 
at its new location warmer than the adjacent solid, so 
in step 4 heat flows from the gas to the solid, cooling 
the gas and causing thermal contraction of the gas to 
a smaller volume. This brings it back to the start of the 
cycle, ready to repeat step 1. 

Although the mass of gas under consideration returns 
to its starting conditions each cycle, its net effect on its 
surroundings is nonzero. First, its thermal expansion oc¬ 
curs at a higher pressure than its thermal contraction, so 
§ p dV > 0: the gas does work on its surroundings, sat¬ 
isfying Rayleigh’s criterion. This work is responsible for 
the positive slope of E versus x in the stack in Fig. 7.2b 
and is represented by the gain element G sta ck^H,i in 
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Fig. 7.2c. All masses of gas within the stack contribute 
to this production of acoustic power; one can think of the 
steady-state operation as due to all masses of gas within 
the stack adding energy to the oscillation every cycle, to 
make up for energy lost from the oscillation elsewhere, 
e.g., in viscous drag in the resonator and acoustic radi¬ 
ation to the surroundings. Second, the gas absorbs heat 
from the solid at the left extreme of its motion, at a rel¬ 
atively warm temperature, and delivers heat to the solid 
farther to the right at a lower temperature. In this way, 
all masses of gas within the stack pass heat along the 
solid, down the temperature gradient from left to right; 
within a single pore, the gas masses are like members of 
a bucket brigade (a line of people fighting a fire by pass¬ 
ing buckets of water from a source of water to the fire 
while passing empty buckets back to the source). This 
transport of heat is responsible for most of H inside the 
stack, shown in Fig. 7.2b. 

This style of engine is called standing wave because 
the time phasing between pressure and motion is close 
to that of a standing wave. (If it were exactly that of 
a standing wave, E would have to be exactly zero at all 
jc.) To achieve the time phasing between pressure and 
volume changes that is necessary for § p dV > 0, im¬ 
perfect thermal contact between the gas and the solid in 
the stack is required, so that the gas can be somewhat 
thermally isolated from the solid during the motion in 
steps 1 and 3 but still exchange significant heat with the 
solid during steps 2 and 4. This imperfect thermal con¬ 
tact occurs because the distance between the gas and 
the nearest solid surface is of the order of <$ t herm> and it 
causes R t herm to be significant, so standing-wave engines 
are inherently inefficient. Nevertheless, standing-wave 
engines are exceptionally simple. They include milli¬ 
watt classroom demonstrations like the illustration in 
Fig. 7.2, similar demonstrations with the addition of 
a water- or air-cooled heat exchanger at the ambient end 
of the stack, research engines [7.25,26] up to several 
kW, and the Taconis and Sondhauss oscillations [7.5,7]. 
Variants based on the same physics of intrinsically irre¬ 
versible heat transfer include the no-stack standing-wave 
engine [7.27], which has two heat exchangers but no 
stack, and the Rijke tube [7.28], which has only a single, 
hot heat exchanger and uses a superposition of steady 
and oscillating flow of air through that heat exchanger 
to create § p dV > 0. 

7.3.2 Traveling-Wave Engines 

One variety of what acousticians call traveling-wave 
engines has been known for almost two centuries as 


a Stirling engine [7.13,29,30], and is illustrated in 
Fig. 7.3a, Fig. 7.3b, Fig. 7.3f, and Fig.7.3g. A regener¬ 
ator bridges the gap between two heat exchangers, one 
at ambient temperature 7 a and the other at hot tempera¬ 
ture 7 h; this assembly lies between two pistons, whose 
oscillations take the gas through a sinusoidal cycle that 
can be approximated as four discrete steps: compres¬ 
sion, displacement rightward toward higher temperature, 
expansion, and displacement leftward toward lower 
temperature. For a small mass of gas in a single pore 
in the heart of the regenerator, the four steps of the cycle 
are illustrated in Fig. 7.3f. In step 1, the gas is com¬ 
pressed by rising pressure, rejecting heat to the nearby 
solid. In step 2, it is moved to the right, toward higher 
temperature, absorbing heat from the solid and experi¬ 
encing thermal expansion as it moves. In step 3, the gas 
is expanded by falling pressure, and absorbs heat from 
the nearby solid. In step 4, the gas is moved leftward, to¬ 
ward lower temperature, rejecting heat to the solid and 
experiencing thermal contraction as it moves. The Stir¬ 
ling engine accomplishes § p dV > 0 in Fig. 7.3g for 
each mass of gas in the regenerator, and this work pro¬ 
duction allows the hot piston to extract more work from 
the gas in each cycle than the ambient piston delivers to 
the gas. 

The similarities and differences between this process 
and the standing-wave process of Fig. 7.2 are instructive. 
Here, the pore size is <^C <$therm, so the thermal contact 
between the gas and the solid in the regenerator is excel¬ 
lent and the gas is always at the temperature of the part 
of the solid to which it is adjacent. Thus, the thermal ex¬ 
pansion and contraction occur during the motion parts 
of the cycle, instead of during the stationary parts of the 
cycle in the standing-wave engine, and the pressure must 
be high during the rightward motion and low during the 
leftward motion to accomplish § p dV > 0. This is the 
time phasing between motion and pressure that occurs 
in a traveling wave [7.14,15], here traveling from left 
to right. The small pore size maintains thermodynami¬ 
cally reversible heat transfer, so Rtherm is negligible, and 
traveling-wave engines have inherently higher efficiency 
than standing-wave engines. One remaining source of 
inefficiency is the viscous resistance R w i sc in the regen¬ 
erator, which can be significant because the small pores 
necessary for thermal efficiency cause R v i sc to be large. 
To minimize the impact of R w i sc , traveling-wave engines 
have |pi| > p m c \U\ | /A, so the magnitude of the spe¬ 
cific acoustic impedance |z a d = \pi \ A/ \ U\ \ is greater 
than that of a traveling wave. 

The gain G listed in Table 7.2 takes on a particularly 
simple form in the tight-pore limit: G reg = AT m /T^ m . 
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Fig. 7.3a-g Traveling-wave engines, (a) A Stirling engine. From left to right , the ambient piston, the ambient heat 
exchanger, the regenerator, the hot heat exchanger, and the hot piston. Time-averaged thermal power Qu is injected 
into the gas at hot temperature 7 h, waste thermal power Qa is removed at ambient temperature 7 a, and net mechanical 
power is extracted by the two pistons, (b) Lumped-element model of the engine in (a), (c) Acoustic-Stirling hybrid 
engine, with the same processes as (a) in the regenerator and its two adjacent heat exchangers, but with additional 
acoustic components replacing the two pistons, (d) Acoustic power E as a function of position x in the device shown 
in (c). Positive power flows in the positive-v direction, so the bottom branch of the curve represents power flowing 
leftward through the feedback path. Total power H is not shown, because it is essentially zero in the regenerator and 
essentially identical to E in the open parts of the device, (e) Lumped-element model of the engine in (c). (f) Close-up 
view of part of one pore in the regenerator of (a) or (c), showing a small mass of gas going through one full cycle, 
imagined as four discrete steps, (g) Plot showing how the pressure p and volume V of that small mass of gas evolve with 
time in a clockwise elliptical trajectory. Tick marks show approximate boundaries between the four steps of the cycle 
shown in (f) 


In the engine of Fig. 7.3a, the initial temperature is 7 a, 
and AT m = Tn~ Ta, so the extra volume velocity that 
the lumped-element model injects at the right end of the 
regenerator is Uaa (7h — 7a) /7a and the total volume 
velocity at the right end of the regenerator is Ua, 1 7h/ 7a . 
Thus, the regenerator acts like an amplifier of volume ve¬ 


locity, with amplification Tr/Ta- If 77 reg is small so that 
pi is nearly the same on both sides of the regenerator, 
then the regenerator amplifies E by nearly 7 h/7a- 
Figure 7.3c-e illustrates a thermoacoustic-Stirling 
hybrid engine [7.16-18], in which the processes in the 
regenerator and its adjacent heat exchangers are the same 
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Fig. 7.4a r b Pulsed combustor, (a) From left to right , check 
valve for admitting fresh air, fuel injection, combustion 
cavity, neck, (b) Lumped-element model of the combustor 


as in the Stirling engine and in Fig. 7.3f and Fig. 7.3g, but 
with acoustic elements replacing the Stirling engine’s 
pistons. The toroidal topology of the thermoacoustic- 
Stirling hybrid allows some of the acoustic power that 
leaves the hot end of the regenerator to be fed back 
to the ambient end of the regenerator, eliminating the 
need for the Stirling engine’s ambient piston. The feed¬ 
back path has an inertance Lf ee d (with an unavoidable 
but small viscous resistance Rf ee d) and a compliance 
Qeed, With 1/LfeedQeed > to 2 and with mL fee d signifi¬ 
cantly smaller than the regenerator viscous resistance 
R reg . These choices let the LfeedQeed feedback path 
boost p\ as acoustic power flows through it, provid¬ 
ing the extra p\ needed to drive U\ into the regenerator 
at its ambient end. 

The Stirling engine in Fig. 7.3a has a hot piston, 

which extracts acoustic power from the gas. From the 

hot piston’s face, mechanical power flows from 7 h to 7 a 

along a temperature gradient in a moving part, either the 

piston itself or a connecting rod. In the thermoacoustic- 

Stirling hybrid engine of Fig. 7.3c, this thermal isolation 

function is accomplished by the thermal buffer tube, 

a thermally stratified column of moving gas. A well- 

designed thermal buffer tube passes acoustic power with 

• • 

little attenuation and has H ~ E, so almost no thermal 

power flows from hot to ambient along the tube and 

little thermal power need be removed at the auxiliary 

heat exchanger at the thermal buffer tube’s ambient end. 

Minimizing streaming and attendant heat convection in 

• • 

thermal buffer tubes, which otherwise causes 77 7 ^ E, is 
a topic of current research in thermoacoustics. 

Traditional Stirling engines are used for propul¬ 
sion in some submarines [7.31] and for auxiliary 
power in boats [7.32], and are under development for 
residential cogeneration of electricity and space heat¬ 
ing [7.33], concentrated solar electricity generation, and 


nuclear generation of electricity for spacecraft [7.34]. 
In most of these applications, piston motion is con¬ 
verted to electricity via relative motion of wires 
and permanent magnets, either with a rotary al¬ 
ternator for crankshaft-coupled pistons or a linear 
alternator (reminiscent of a loudspeaker) for resonant 
free-piston configurations. Thermoacoustic-Stirling hy¬ 
brid engines are under consideration for small-scale 
natural-gas liquefaction [7.25, 35] and for spacecraft 
power [7.36]. In the former case, acoustic power 
is fed directly from the engine to cryogenic acous¬ 
tic refrigerators, without transduction to electrical 
power. 

7.3.3 Combustion 

In the standing-wave and traveling-wave engines, 
Rayleigh’s criterion § p dV > 0 is met with volume 
changes that arise from temperature changes; those 
temperature changes, in turn, arise from thermal con¬ 
tact between the gas and nearby solid surfaces. In 
pulsed combustion, the volume changes needed to meet 
Rayleigh’s criterion arise from both temperature and 
mole-number changes, which in turn are due to time- 
dependent chemical reactions whose rate is controlled 
by the time-dependent pressure or time-dependent ve¬ 
locity [7.37,38]. 

Figure 7.4 illustrates one configuration in which 
pulsed combustion can occur. At the closed end of 
a closed-open resonator, a check valve periodically lets 
fresh air into the resonator and a fuel injector adds fuel, 
either steadily or periodically. If the rate of the exother¬ 
mic chemical reaction increases with pressure (e.g., via 
the temperature’s adiabatic dependence on pressure), 
positive dV occurs when p is high, meeting Rayleigh’s 
criterion. A four-step diagram of the process, resem¬ 
bling Fig. 7. 2d and Fig. 7.3f, is not included in Fig. 7.4 
because the process is fundamentally not cyclic: a given 
mass of gas does not return to its starting conditions, 
but rather each mass of fresh-air-fuel mixture burns and 
expands only once. 

Combustion instabilities can occur in rockets, jet 
engines, and gas turbines, with potentially devastat¬ 
ing consequences if the pressure oscillations are high 
enough to cause structural damage. Much of the 
literature on thermoacoustic combustion is devoted to 
understanding such oscillations and using active or pas¬ 
sive means to prevent them. However, some devices such 
as high-efficiency residential gas-fired furnaces deliber¬ 
ately use pulsed combustion as illustrated in Fig. 7.4 to 
pump fresh air in and exhaust gases out of the combus- 
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tor. This eliminates the need to leave the exhaust gases fraction of the heat of combustion can be delivered to 
hot enough for strong chimney convection, so a larger the home. 


7.4 Dissipation 


The dissipative processes represented above by R w i sc 
and ^ therm occur whenever gas-borne sound inter¬ 
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acts with solid surfaces. Figure 7.5 illustrates this 
in the case of a short length dx of a large-radius 
duct with no axial temperature gradient. The ori¬ 
gin of the viscous dissipation of acoustic power is 
viscous shear within the viscous penetration depth 
<$vi sc , as shown in Fig. 7.5c. Most people find vis¬ 
cous dissipation intuitively plausible, imagining the 
frictional dissipation of mechanical energy when 
one surface rubs on another. More subtle is the 
thermal relaxation of acoustic power, illustrated in 
Fig. 7.5d and Fig. 7.5e. Gas is pressurized nearly 
adiabatically in step 1, then shrinks during ther¬ 
mal equilibration with the surface in step 2. It 
is depressurized nearly adiabatically in step 3, and 
then thermally expands during thermal equilibra¬ 
tion with the surface during step 4. As shown in 
Fig. 7.5e, the net effect is § p dV < 0: the gas ab¬ 
sorbs acoustic power from the wave, because the 
contraction occurs at high pressure and the ex¬ 
pansion at low pressure. To avoid a hopelessly 
cluttered illustration, Fig.7.5d shows the thermal- 
hysteresis process superimposed on the left-right 
oscillating motion in steps 1 and 3, but the thermal- 
hysteresis process occurs even in the absence of such 
motion. 

Differentiating (7.18) with respect to x shows that 
the dissipation of acoustic power in the duct in length 

Fig. 7.5a-e Boundary dissipation in acoustic waves, 
(a) A duct with no temperature gradient, with one short 
length dx identified, (b) Each length dx has inertance, 
viscous resistance, compliance, and thermal hysteresis re¬ 
sistance. (c) The dissipation of acoustic power by viscous 
resistance is due to shear in the gas within roughly 5 v i sc 
of the boundary, here occurring during steps 1 and 3 of 
the cycle, (d) and (e) The dissipation of acoustic power 
by thermal relaxation hysteresis occurs within roughly 
^therm of the boundary. Gas is pressurized nearly adiabat¬ 
ically in step 1, then shrinks during thermal equilibration 
with the surface in step 2. It is depressurized nearly adi¬ 
abatically in step 3, and then thermally expands during 
thermal equilibration with the surface during step 4. The 
net effect is that the gas absorbs acoustic power from 
the wave 
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dx is given by 

1 (dpi dU\ \ 

d£=-Re(-^-t/i+pi- r L dx , ( 7 . 21 ) 

2 \ dx dx J 

and examination of (7.10) and (7.9) associates R w i sc 
with the first term and R t herm with the second term. 
Expressions for 7? v i sc and 7? t herm in the boundary-layer 
approximation are given in Table 7.2, and allow the ex¬ 
pression of (7.21) in terms of the dissipation of acoustic 

7.5 Refrigeration 

7.5.1 Standing-Wave Refrigeration 

The thermal-hysteresis process described in Sect. 7.4 
consumes acoustic power without doing anything 
thermodynamically useful. Standing-wave refrigeration 
consumes acoustic power via a similar process, but 
achieves a thermodynamic purpose: pumping heat up 
a temperature gradient, either to remove heat from 
a temperature below ambient or (less commonly) to de¬ 
liver heat to a temperature above ambient. Figure 7.6a-c 
shows a standing-wave refrigerator of the style pi¬ 
oneered by Hofler [7.8, 9] and recently studied by 
Tijani [7.39]. At the left end, a driver such as a loud¬ 
speaker injects acoustic power E , which flows rightward 
through the stack, causing a leftward flow of total energy 
77 . 

The most important process in a standing-wave re¬ 
frigerator is illustrated in Fig. 7.6d and Fig. 7.6e from 
a Fagrangian perspective. Fig. 7.6d shows a greatly mag¬ 
nified view of a small mass of gas inside one pore of 
the stack. The sinusoidal thermodynamic cycle of that 
mass of gas in pressure p and volume V is shown in 
Fig. 7.6e; the mass’s temperature, entropy, density, and 
other properties also vary sinusoidally in time. How¬ 
ever, for qualitative understanding of the processes, we 
describe them as if they are a time series of four discrete 
steps, numbered 1-4 in Fig. 7.6d and Fig. 7.6e. In step 1, 
the gas is simultaneously compressed to a smaller vol¬ 
ume and moved leftward by the wave. Thermal contact 
is imperfect in the pores of a stack, so during step 1 the 
gas experiences a nearly adiabatic temperature increase 
due to the pressure increase that causes the compression. 
It arrives at its new location warmer than the adjacent 
solid because the temperature gradient in the solid is less 
than the critical temperature gradient Vf cr i t defined in 


power per unit of surface area S : 

• 1 2 

d E =—Pm \u\ | <u<5yj sc dS 

1 I 12 

+ T (V ~ 1) ~~2 ^^therm dS , (7.22) 

4 p m c z 

where u\ = U\/ A is the velocity outside the boundary 
layer, parallel to the surface. Each term in this result ex¬ 
presses a dissipation as the product of a stored energy 
per unit volume p m |^il 2 /4 or (y — 1 ) \p\ | 2 /4p m c 2 , 
a volume <$ v i sc dS or <$ t herm d S, and a rate 00 . 


(7.20). Thus, in step 2, heat flows from the gas into the 
solid, cooling the gas and causing thermal contraction 
of the gas to a smaller volume. In step 3, the gas is si¬ 
multaneously expanded to a larger volume and moved 
rightward by the wave. It arrives at its new location 
cooler than the adjacent solid, so in step 4 heat flows 
from the solid to the gas, warming the gas and causing 
thermal expansion of the gas to a larger volume. This 
brings it back to the start of the cycle, ready to repeat 
step 1 . 

The mass of gas shown in Fig. 7.6d has two time- 
averaged effects on its surroundings. First, the gas 
absorbs heat from the solid at the right extreme of its mo¬ 
tion, at a relatively cool temperature, and delivers heat 
to the solid farther to the left at a higher temperature. 
In this way, all masses of gas within the stack pass heat 
along the solid, up the temperature gradient from right to 
left—within a single pore, the gas masses are like mem¬ 
bers of a bucket brigade passing water. This provides the 
desired refrigeration or heat-pumping effect. If the left 
heat exchanger is held at ambient temperature, as shown 
in Fig. 7.6a, then the system is a refrigerator, absorb- 
ing thermal power Qc from an external load at the right 
heat exchanger at one end of the bucket brigade at 7c, 
as waste thermal power is rejected to an external ambi¬ 
ent heat sink at the other end of the bucket brigade at 7 a. 
(The system functions as a heat pump if the right heat ex¬ 
changer is held at ambient temperature; then the left heat 
exchanger is above ambient temperature.) Second, the 
gas’s thermal expansion occurs at a lower pressure than 
its thermal contraction, so (f p dV < 0: the gas absorbs 
work from its surroundings. This work is responsible for 
the negative slope of E versus x in the stack in Fig. 7.6b 
and is represented by the gain element G sta ck^A,i in 
Fig. 7.6c. All masses of gas within the stack contribute 
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to this consumption of acoustic power, which must be 
supplied by the driver. 

As in the standing-wave engine, the time phasing be¬ 
tween pressure and motion in the standing-wave refriger¬ 
ator is close to that of a standing wave. Imperfect thermal 
contact between the gas and the solid in the stack is re- 
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quired to keep the gas rather isolated from the solid dur¬ 
ing the motion in steps 1 and 3 but still able to exchange 
significant heat with the solid during steps 2 and 4. This 
imperfect thermal contact occurs because the distance 
between the gas and the nearest solid surface is of the 
order of <$ t herm, and it causes 7? t herm to be significant, so 
standing-wave refrigerators are inherently inefficient. 

7.5.2 Traveling-Wave Refrigeration 

Several varieties of traveling-wave refrigerator are com¬ 
mercially available or under development. At their 
core is a regenerator, in which the process shown in 
Fig. 7.7a, b operates. In step 1 of the process, the gas 
is compressed by rising pressure, rejecting heat to the 
nearby solid. In step 2, it is moved to the right, toward 
lower temperature, rejecting heat to the solid and ex¬ 
periencing thermal contraction as it moves. In step 3, 
the gas is expanded by falling pressure, and absorbs 
heat from the nearby solid. In step 4, the gas is moved 
leftward, toward higher temperature, absorbing heat 
from the solid and experiencing thermal expansion as 
it moves. The heat transfers between gas and solid in 
steps 2 and 4 are equal and opposite, so the net ther¬ 
mal effect of each mass of gas on the solid is due to 
steps 1 and 3, and is to move heat from right to left, 
up the temperature gradient. As before, the motion of 
any particular mass of gas is less than the length of 
the regenerator, so the heat is passed bucket-brigade 
fashion from the cold end of the regenerator to the am¬ 
bient end. Each mass of gas absorbs § p dV of acoustic 
power from the wave as shown in Fig. 7.7b, because 
the thermal contraction in step 2 occurs at high pres- 

Fig. 7.6a-e A standing-wave refrigerator, (a) From left to 
right, a piston, ambient heat exchanger, stack, cold heat 
exchanger, tube, and tank. Acoustic power is supplied to 
the gas by the piston to maintain the standing wave, and 
results in thermal power Qc being absorbed by the gas from 
a load at cold temperature Tq while waste thermal power 
Qa is rejected by the gas to an external heat sink at ambient 
temperature 7 a- (b) Total power flow H and acoustic power 
E as functions of position x in the device. Positive power 
flows in the positive-v direction, (c) Lumped-element model 
of the device. (d) Close-up view of part of one pore in the 
stack of (a), showing a small mass of gas going through one 
full cycle, imagined as four discrete steps, (e) Plot showing 
how the pressure p and volume V of that small mass of gas 
evolve with time in a counter-clockwise elliptical trajectory. 
Tick marks show approximate boundaries between the four 
steps of the cycle shown in (d) 
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Fig.7.7a-g Traveling-wave refrigerators, (a) Close-up view of part of one pore in the regenerator of a traveling-wave 
refrigerator, showing a small mass of gas going through one full cycle, imagined as four discrete steps, (b) Plot showing 
how the pressure p and volume V of that small mass of gas evolve with time in a counterclockwise elliptical trajectory. 
Tick marks show approximate boundaries between the four steps of the cycle shown in (a), (c) A traditional Stirling 
refrigerator. From left to right , the ambient driver piston; the ambient heat exchanger, the regenerator, and the cold heat 
exchanger are across the top; the displacer is below. Time-averaged thermal power Qc is removed by the gas from an 
external heat load at cold temperature 7c, while waste thermal power Qa is rejected by the gas to an external heat sink at 
ambient temperature 7\ and net acoustic power is supplied to the gas by the driver, (d) Pulse-tube refrigerator, with the 
same processes as (c) in the regenerator and its two adjacent heat exchangers, but with the displacer replaced by acoustic 
components, (e) Acoustic power E and total power H as functions of position x in the device shown in (d). Positive 
power flows in the positive-v direction, (f) Lumped-element model of the pulse-tube refrigerator in (d). (g) A refrigerator 
in which the thermal-buffer function and moving-mass function of (c)’s displacer are in two separate components 


sure and the thermal expansion in step 4 occurs at low 
pressure. The small pore size, <^C <$therm, maintains 
thermodynamically reversible heat transfer, so 7? t herm is 
negligible, and traveling-wave refrigerators have an in¬ 
herently high efficiency. Acoustically, the process shown 
in Fig. 7.7a represents acoustic power traveling from left 


to right through the regenerator, and being partly con¬ 
sumed as it goes. Different varieties of traveling-wave 
refrigerator use different methods to create the neces¬ 
sary amplitudes and relative time phasing of motion and 
pressure to achieve this acoustic power flow through the 
regenerator. 
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The traditional Stirling refrigerator [7.40] uses two 
moving pistons, which can be either crankshaft-coupled 
in a configuration like that of the engine in Fig. 7.3a or of 
th efree-displacer variety shown in Fig. 7.7c. In the free- 
displacer Stirling refrigerator (often called free piston ), 
the solid displacer moves in response to the gas-pressure 
forces on it, without linkage to any external motor. Its 
area and mass are selected to give its motion the de¬ 
sired amplitude and time phase [7.29]. Acoustic power 
is transmitted through it, from right to left, so that the 
acoustic power that flows out of the right end of the re¬ 
generator is fed through the displacer to the left, added 
to the acoustic power supplied by the driver, and injected 
into the left end of the regenerator. The displacer must 
be a thermal insulator, because its right end is at Tq and 
its left end at 7 a. Crankshaft-coupled Stirling refriger¬ 
ators were used in the 19th century to keep beef cold 
on the long sea voyage from South America to Britain. 
Free-piston Stirling cryocoolers are in common use to¬ 
day for cooling infrared sensors in military night-vision 
goggles and surveillance satellites, and a free-piston Stir¬ 
ling refrigerator built into a small, portable ice chest is 
commercially available at low cost for picnics [7.41]. 

The pulse-tube refrigerator [7.12], illustrated in 
Fig.7.7d-f, uses only one piston. The acoustic power 
flowing out of the right end of the regenerator is absorbed 
in an acoustic impedance, instead of being fed back to 
the left end of the regenerator. This gives the pulse-tube 


7.6 Mixture Separation 

In thermoacoustic mixture separation, acoustic power 
causes the components of a gas mixture to separate [7.47, 
48]. The process is loosely analogous to the pumping of 
heat through the stack in a standing-wave refrigerator. 
The expenditure of acoustic power results in an increase 
in the Gibbs free energy of the mixture’s components, 
and the efficiency of the process is comparable to that 
of some other practical separation processes [7.49]. 

Figure 7.8 illustrates the process for a binary gas 
mixture of heavy and light molecules, whose motions 
are indicated by filled and open arrows, respectively. 
The motion of the molecules in steps 1 and 3 of the 
process is bulk motion of the gas, and the motion 
during steps 2 and 4 is thermal diffusion, in which 
light molecules diffuse toward higher temperature and 
heavy molecules diffuse toward lower temperature. In 
step 1, the gas in the center of the pore moves left¬ 
ward and its pressure rises, so its temperature rises 


refrigerator a lower efficiency than the Stirling refrigera¬ 
tor, but for cryogenic applications the reduced efficiency 
is a small price to pay for the elimination of the Stir¬ 
ling’s cold moving part. The lumped-element model of 
the impedance and adjacent tank shown in Fig. 7.7f cor¬ 
rectly suggests that proper design of L impe d, limped, 
and Qank can create almost any desired ratio of | U\ | to 
\p\ \ and time phasing between U\ and p\ at the right 
end of the regenerator. The desired impedance is often 
achieved by choosing ^ visc ~ coL[ m pe d » l/o>Ci mpe d. 
The so-called pulse tube, and the ambient heat ex¬ 
changer to its right, thermally isolate the cold heat 
exchanger from the dissipation of acoustic power in 
the impedance, just as the thermal buffer tube performs 
that function in Fig. 7.3c. Pulse-tube refrigerators are 
in common use today in satellites [7.42] and are under 
development for many other applications such as small- 
scale oxygen liquefaction and cooling superconducting 
equipment. 

Another variation of the theme is shown in Fig. 7.7g, 
where the thermal buffer tube and the moving mass 
perform the thermal-insulation and inertial functions 
of the free displacer of Fig. 7.7c in two separate lo¬ 
cations [7.43,44]. This variety is under development for 
commercial food refrigeration [7.45]. Yet another vari¬ 
ety [7.46] replaces the inertial moving mass of Fig. 7.7g 
with inertial moving gas, similar to the thermoacoustic- 
Stirling hybrid engine of Fig. 7.3c. 


nearly adiabatically. Trapped by the viscous bound¬ 
ary layer near the wall of the tube is some other gas 
that does not move and whose thermal contact with 
the pore wall keeps its temperature from rising. During 
step 2, the temperature difference between the central 
gas and the peripheral gas causes light atoms to dif¬ 
fuse into the center and heavy atoms to diffuse out of 
the center, enriching the center in light atoms. Bulk 
motion rightward in step 3 carries this central light- 
enriched gas to the right. Simultaneously, the pressure 
drops in step 3, and hence the central temperature drops 
nearly adiabatically. Thus thermal diffusion in step 4 
pulls heavy molecules into the center and drives light 
molecules out of the center, leaving the center enriched 
in heavy molecules so that the leftward motion dur¬ 
ing step 1 carries heavy-enriched gas leftward. The net 
effect of steps 1 and 3 is to move heavy molecules 
leftward and light molecules rightward. The process 
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takes place when the mole-fraction gradient is less than 
a critical gradient, analogous to the critical tempera¬ 
ture gradient below which standing-wave refrigeration 
occurs. 

As is evident from Fig. 7.8, the process works best 
for tubes with r h somewhat larger than <$ t herm and <S v i sc . 
However, unlike the description above, the process ac- 


Fig.7.8 Close-up view of part of one pore in a thermo¬ 
acoustic mixture separator, showing a small body of gas in 
the center of the pore going through one full cycle, imagined 
as four discrete steps, and exchanging mass with neighbor¬ 
ing, immobile gas. The gray and brown arrows in steps 2 
and 4 signify thermally driven diffusion of light and heavy 
molecules between the central gas and the peripheral gas 
close to the wall of the pore. In steps 1 and 3, the gray 
and brown arrows signify bulk motion of the gas, without 
diffusion 

tually works best near traveling-wave phasing, because 
two 45° phase shifts were ignored in that description. 
High separation purities require long tubes, and high 
mole fluxes will require many passages in parallel, per¬ 
haps in structures similar to the stacks of standing-wave 
engines and refrigerators. 

In a 2 m-long tube, a 50-50 helium-argon mixture 
has been separated to yield 70% helium, 30% argon at 
one end and 30% helium, 70% argon at the other end; 
and a measurable enrichment of 22 Ne from natural neon 
has been achieved [7.50]. 
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